Photon transfer in ultrastrongly coupled three-cavity arrays 



CO 

o 

CN 

u 

Oh 
< 

en 

CN 



Oh 
i 

■4— > 

0< 



CN 
CN 

O 
en 

> 

X 



S. Felicetti, 1 G. Romero, 1 D. Rossini, 2 R. Fazio, 2,3 and E. Solano 1 ' 4 

1 Departamento de Quimica Fisica, Universidad del Pais Vasco UPV/EHU, Apartado 644, 48080 Bilbao, Spain 

2 NEST, Scuola Norrnale Superiore and Istituto di Nanoscienze - CNR, Pisa, Italy 

Centre for Quantum Technologies, National University of Singapore, Republic of Singapore 

IKERBASQUE, Basque Foundation for Science, Alameda Urquijo 36, 48011 Bilbao, Spain 

(Dated: April 24, 2013) 

We study the photon transfer along a linear array of three coupled cavities where the central one contains an 
interacting two-level system in the strong and ultrastrong coupling regimes. We find that an inhomogeneously 
coupled array forbids a complete single-photon transfer between the external cavities when the central one 
performs a Jaynes-Cummings dynamics. This is not the case in the ultrastrong coupling regime, where the 
system exhibits singularities in the photon transfer time as a function of the cavity-qubit coupling strength. 
Our model can be implemented within the state-of-the-art circuit quantum electrodynamics technology and it 
represents a building block for studying photon state transfer through scalable cavity arrays. 



Introduction. — Light-matter interaction controls some of 
the most fundamental processes in nature QJ. Its study has 
allowed the development of impressive architectures where 
fundamentals of quantum mechanics can be tested |J2). The 
high level of control achieved on these setups [3-5], moti- 
vated an increasing interest in the study of strongly correlated 
systems and collective phenomena. In particular, photon lat- 
tice models involving coupled cavity arrays have been devised 
in such a way that each cavity interacts with a single two-level 
system (qubit), thus realizing the Jaynes-Cummings-Hubbard 
model [6-8]. The possibility of implementing different ge- 
ometries enables one to engineer quantum networks for dis- 
tributed quantum information processing [9]. These lattice 
models have proved useful to describe the scattering of a 
single-photon interacting with a qubit in a one-dimensional 
waveguide [10|. A similar approach can be used to con- 
sider homogeneously coupled cavity arrays [11], where the 
central cavity-qubit interaction is in the strong coupling (SC) 
regime. However, when the rotating-wave approximation 
(RWA) ceases to be valid, the single-photon scattering prop- 
erties have been only recently studied [12]. Note that cir- 
cuit QED technologies lfT3lfT5l . where photon lattice models 
have been proposed |[T6l - n"8"l . include disorder in the cavity 
coupling distribution due to imperfections. In addition, cir- 
cuit QED has allowed the advent of the ultrastrong coupling 
(USC) regime of light-matter interaction |fT9H2TI . described 
by the quantum Rabi model 11221 . The latter displayed impor- 
tant consequences in strongly correlated systems, where a Z2 
parity-breaking quantum phase transition was predicted [23|. 

In this Letter, we consider the problem of photon trans- 
fer in a linear array of three coupled cavities, where a two- 
level system interacts at the central site in the SC and USC 
regimes (see Fig. [T}. This configuration can be thought as a 
microwave analogue |24| of the superconducting Josephson 
interferometer. We also include disorder in the cavity-cavity 
coupling, thus mimicking experimental imperfections. Un- 
der these conditions, we are able to unveil the following fea- 
tures: (i) in the SC regime, and for finite values of the hop- 
ping amplitudes and of the cavity-qubit coupling, a single ex- 
citation initially localized in the leftmost cavity is strictly for- 



bidden to fully populate the rightmost cavity, similar to the 
delocalization-localization phenomenon 1251 : (ii) in the USC 
regime, the above restriction does not hold any more, and a 
single-excitation can fully populated the rightmost cavity for 
almost arbitrary Hamiltonian parameters; (iii) if the qubit tran- 
sition energy is null, we find a complete and periodic excita- 
tion transfer from the leftmost to the rightmost cavity, inde- 
pendent of the cavity-qubit coupling strength. Our scheme 
represents a feasible building block |fT6T[T8l [26 1 to study pho- 
ton excitation and state transfer towards scalable cavity arrays. 
The model. — Our model consists of an array of three single- 
mode cavities, where the central site interacts with a two-level 
system in the SC or in the USC regime. A schematic repre- 
sentation is shown in Fig.fTlwhere each cavity is linked with 
its neighbor through a hopping interaction that, in general, is 
not weak enough to consider the RWA. The corresponding 
Hamiltonian reads 



H = ^wea\a e + -^-a z + ga x (al+a 2 ) 



(=1 
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^ Jn{a\ai + i + a\a\ +1 + H.c), 
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where ag(a e ) is the annihilation (creation) operator for pho- 
tons on the f-th cavity (£ = 1,2, 3), u>i being the characteristic 




FIG. 1. (Color online). Linear chain of three microwave cavities. 
The central cavity is coupled to a two-level quantum system in the 
strong coupling or ultrastrong coupling regime. 




FIG. 2. (Color online). Oscillation amplitude T defined in Eq. pi 
as function of hopping parameter J2 and for different parameter sets. 
In panel (a), for cavity-qubit coupling g = 0, we display T for J\ = 
0.001 (continuous blue line), J\ — 0.005 (dashed black line), and 
Ji = 0.01 (dot-dashed red line). In panel (b), for J x = 0.001, we 
consider g — (continuous blue line), g = 0.002 (dashed black 
line), and g — 0.01 (dot-dashed red line). 



frequencies and Ji are nearest-neighbor hopping amplitudes. 
The qubit of frequency u q is located inside the central cavity 
and is described by the Pauli matrices a a (a — x, y, z), while 
g denotes the cavity-qubit coupling strength. The complexity 
of the Hamiltonian in Eq. (fTh is associated with the appearance 
of counter-rotating terms in the cavity-qubit and cavity-cavity 
interaction. Hereafter, we consider identical cavities (ug = u), 
and the condition Jp/u ~ 0.1 to assure a faithful numerical 
analysis of Hamiltonian (fTh. We also set the energy scales in 
units of the resonator frequency ui. 

The system is initialized in the state |V>o) = 1 100) <8> |g), 
corresponding to having a single photon in the leftmost cav- 
ity, zero in the others, and the qubit in its ground state. We 
address the dynamics dictated by the Hamiltonian (jl), and 
study the single-excitation transfer along our three-cavity ar- 
ray. Depending on the ratio g/u), two regimes can be identi- 
fied: the SC regime for g/u < 0.1, and the USC regime for 
0.1<p/w<l. 

Single-photon transfer in SC regime. — When g/u <C 1, the 
RWA provides a faithful description of the cavity-qubit dy- 
namics, so that we can neglect the counter-rotating terms in 
Eq. i\\: a x {a\ + 02) — > (o~+a,2 + a^a 2 )- In the regime 
where hoppings Ji/u <C 1, the cavity-cavity RWA also holds, 
and the U(\) symmetry provides the conservation of the to- 
tal number of excitations. In this case, the time-evolution 
of the system will necessarily lead to a state of the form 
|V(t))=a J000) ® |e) + (/3 |100) +7 |010) +8 1001) ) <8> |g). A 
full analytical solution can be found in the interaction picture 
and directly solving the Schrodinger equation. At resonance, 
where qubit and resonator frequencies coincide (u q = u), we 



can derive explicitly the probability amplitude for finding a 
photon in the rightmost cavity: S(t) = T [cos (At) — 1] /2, 
where A = \Jg 2 + J 2 + J\ and the amplitude reads 

2JiJ 2 



T 



Jl 



(2) 



9 1 + Jl 

The above result unveils a competition between cavity-cavity 
and cavity-qubit interaction. Figure [2k indeed shows that, in- 
troducing disorder in the cavity-cavity couplings ( J\ 7^ J2), 
the single-photon transfer has a counter-intuitive dependence 
on the hopping terms. In particular, for given parameters 
g, u, and J\, the excitation transfer to the rightmost cavity 
exhibits a nonmonotonic behavior with increasing J2 and is 
maximum for J\ = J2. Note that only in the homogeneous 
case (Ji = J2) and for a negligible cavity-qubit coupling g, 
it is possible to have photon transfer with unit probability. On 
the other hand, for any finite value of g, the latter is strictly 
forbidden. We also notice that, in Fig. Eh and for fixed J\, 
increasing g decreases the amplitude T. This behavior has 
a simple explanation: when a single excitation tries to move 
from the left to the central cavity, it is scattered back by the 
cavity-qubit system without being fully absorbed. In fact, the 
probability for exciting the qubit \a\ 2 is inversely proportional 
to the square of the coupling strength. 

Single-photon transfer in the USC regime. — If the coupling 
strength g and the resonator frequency u satisfy 0.1 <g/u> < 
1, the system enters the USC regime. In this case, photons 
are spontaneously generated from the vacuum such that the 
total number of excitations grows with the ratio g/u, enlarg- 
ing unavoidably the associated Hilbert space 11271 . In order 
to provide a reliable system dynamics, we performed a fourth 
order- Trotter expansion of the evolution operator [28 , 29 1. 

Figure [3] shows the time-evolution of the average photon 
number in each cavity, N^(t) = (ip(t)\a]ai\?p(t)), starting 
from the state |^o) that evolves according to Hamiltonian (fTJi. 
At first glance, one recognizes a highly irregular behavior of 
A^pL (t), which arises from the counter-rotating terms in the 
cavity-qubit interaction. Remarkably, this is developed by the 
unitary evolution of the system itself, and is not due to the 
limited time-resolution of our simulations. In order to quan- 
tify this behavior, we analyse the fractal dimension of N^ h (t) 
by using the modified box counting algorithm |30|. This con- 
sists in dividing the total time interval in segments of size r, 
and then covering the data with a set of rectangular boxes 
of size tx Aj (Aj is the largest excursion of the curve in 
the i-th region r). One then computes the average excursion 
M(t) — J2i Aj/t. The dimension T> of the curve is defined 
by T> = — log T M(t). One finds T> — 1 for a straight line, 
and T> = 2 for a periodic curve. Indeed, for times much larger 
than the period, a periodic curve uniformly covers a rectangu- 
lar region. Any value of T> between these integer values en- 
tails the fractality of the curve [31 1. In the case of Fig. [3^, we 
obtained T> w 1.66 (inset to the figure). Furthermore, increas- 
ing g, we found a fractal dimension which rapidly decreases 
from T> ~ 2 (quasi-periodic curve) down to non integer values 
close to T> ~ 1.5 for large cavity-qubit couplings. 
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FIG. 3. (Color online). Average photon number in each cav- 
ity AK;(t) at resonance condition (Uq — u), and for homogeneous 
cavity-cavity couplings (Ji = J2=0.1). The red line stands for the 
leftmost cavity, the blue line for the rightmost cavity, while the green 
line for the central cavity. The main panels refer to different cavity- 
qubit couplings: g = 0.9 (a), g = 0.85 (b). The inset shows the 
box counting analysis for jV2 h (t) (g — 0.9), and displays M as a 
function of r (see the main text and Ref. 1311 ). A power-law fit of 
the intermediate region gives a fractal dimension T> ~ 1.66. 



Two fundamental differences with respect to the SC regime 
can be found in the excitation transfer analysis. First, in the 
USC regime a single photon excitation can be completely 
transferred from the leftmost to the rigthmost cavity, also for 
a finite value of the coupling strength g [see, e.g., Fig. 13k], 
This is the typical situation for almost arbitrary Hamiltonian 
parameters in the USC regime, a disparity that can be ex- 
plained as follows. In the SC regime, the initial state \ipo) 
has a finite overlap with the eigenstate \E) — (g\ 100) <8> |g) + 
J|000) ® |e))/\/g 2 + J 2 , which is conserved by the time- 
evolution. Consequently, the state |001) <8> |g) is not accessi- 
ble (as far as g ^ 0). On the other hand, in the USC regime 
the state \E) is no longer a system eigenstate. In fact, any 
eigenstate of Hamiltonian (fT) has a structure that covers the 
full system Hilbert space, thus enhancing the overlap with the 
state 1 001) ® |g). Consequently, complete transfer to the right- 
most cavity is generally allowed. 

The second feature to be highlighted is that the photon 
transfer is strongly inhibited for a specific value of the cavity - 
qubit coupling strength, as displayed in Fig. [3J5. In particular, 
in Fig. HI we analyzed the inversion time T; nv , defined as the 
time where the average photon number becomes bigger in the 
rightmost cavity than in the leftmost cavity, as a function of 
the coupling strength g, and for different values of the cavity- 
cavity coupling strength. We found that there exists a critical 
value g c , dependent on the system parameters, for which the 
time needed to observe single-photon transfer dramatically in- 
creases. For our set of parameters, fixing J\ =0.1 we found 
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FIG. 4. (Color online). Population inversion time Ti nv as a function 
of the cavity-qubit coupling constant, defined as the time in which 
Wph (*) becomes bigger than N$(t) [that is, A/V ph (T inv ) = in 
Eq. fib]. Here, we have fixed J\ = 0.1. Note the SC-USC transition 
for g « 0.14 and the inhibition of state transfer for g = g c « 0.85. 
The inset displays the maximum value of the inversion time T ir ^ 
that is reached at g c , as a function of J2 1291 . 



g c « 0.94 — 0.97 x J2, while the inversion time T^™ in cor- 
respondence to such value exhibits a quite irregular pattern of 
oscillations with J 2 (inset to Fig. El. To our knowledge, this is 
the first observation of this behavior in the excitation transfer 
properties of a cavity QED system. This phenomenon is spe- 
cific to the USC dynamics and occurs in the higher-coupling 
region of the Rabi model (g/ui > 0.4) ll32l . a zone where the 
photon production exceeds the RWA predictions, and an ana- 
lytical treatment becomes difficult despite the integrability of 
the model [22|. 

We highlight the sudden increase of T- mv occurring in the 
lower-coupling region of the Rabi model, which is zoomed-in 
in Fig. 5. This abrupt behavior is due to the SC/USC regime 
crossover. If the RWA holds, the population inversion time 
is given by: T inv = arccos(l — A 2 /(J 2 + J|)) /A, where A 
has been defined previously. For g > \J J 2 + J 2 , the popula- 
tion inversion never occurs. Contrariwise, when the counter- 
rotating terms are taken into account, we observe population 
inversion at finite time also for larger value of the coupling g 
For our three-cavity array, the parameter g t = g/s/Jf 
provides us with an operational definition for the SC/USC 
transition of the cavity-qubit interaction. For small values of 
gt, the RWA describes correctly the transfer dynamics. When 
gt > 1, the full model behavior differs quantitatively and qual- 
itatively from the JC model predictions. We observe also that 
varying J 2 with respect to J\ results in a longer population in- 
version time: in both regimes, inhomogeneity in the hopping 
terms hinders the excitation transfer. 

We have also analyzed the transfer of coherent states, 
\§oi=\a), and of arbitrary linear superposition states |</> ) = 
p |0) + e l6 \J\ — p 2 1 1) , with p randomly chosen in the interval 
[0,1]. We numerically simulated the system evolution, setting 
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FIG. 5. (Color online). Population inversion time T\ nv as a function 
of the cavity-qubit coupling constant g, and for fixed value 3\ = 0.1. 
The blue and the black lines correspond to the homogeneous case. 
The blue line corresponds to the analytical solution when the RWA 
holds for cavity-cavity and cavity-qubit interaction. The black, red 
and green lines are obtained through numerical simulations of the 
full model ruled by the Hamiltonian (T). 



IV'o) = |0o ) |0) |0) |g) as initial state, and we recorded the 
behavior of the transfer fidelity, defined as: F — Tr[pop(t)], 
where p — |0 O ) (4>o\ an d P{p) i s me state of the rightmost 
cavity at time t. The results we found are consistent with those 
relative to the case of single-excitation transfer. The interplay 
between the hopping constant J and the cavity-qubit coupling 
strength g, rules the state transfer dynamics: with increasing 
J the transfer is more likely to happen, while increasing g re- 
sults in smaller values of the transfer fidelity. Specifically, in 
the SC regime, as far as g ^ 0, it is impossible to observe 
complete state transfer, i.e. F < 1 at any time. This is not 
the case in the USC regime, where the transfer fidelity, be 
in a linear superposition or a coherent state, can be close to 
unity, also when g and J are of the same order. These results 
show the peculiar features of state transfer physics beyond the 
RWA, and may pave the way for developing a general theory 
in presence of USC regimes. Further details on state transfer 
features can be found in the supplemental material. 

Degenerate qubit case. — When the qubit frequency van- 
ishes, a closed analytical solution of the dynamics governed 
by Hamiltonian (fill is available if we consider the RWA in 
the cavity-cavity interaction. We may consider the quantum 
simulation of this model by means of a coupled cavity-qubit 
system in the SC regime, see Fig. [Tland the application of a 
strong classical driving to the qubit l33l[34l . It can be shown 
that the dynamics will be ruled by the effective Hamiltonian 
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(3) 

In this case, the excitation transfer exhibits a smooth periodic 
behavior in time, which is independent of the cavity-qubit 
coupling strength, and allows a complete transfer at regular 



times. In fact, the time evolution of the difference between 
the average photon number in the leftmost and rightmost cav- 
ities, starting from the initial state \ipo) = 1 100) (8> |g) in the 
homogeneous case, reads 



(3), 



AiV ph (t) = N$>{t) - Ng>(t) = cos (V2Jt) 



(4) 



This result holds both in the SC and in the USC regime, 
hence the counter-rotating terms do not modify the excitation 
transfer properties of the system. We point out that, when 
g/u>> 0.1, spontaneous photon generation occurs: (a[ai) and 
(0303) keep a fractal time dependence, despite their difference 
follows a smooth and regular behavior. 

Conclusions and outlook. — We studied excitation transfer 
in an array of three coupled cavities, where a two-level sys- 
tem interacts with the central cavity, focusing on the com- 
parison between strong and ultrastrong coupling regimes of 
light-matter interaction. In the SC regime, the cavity-qubit 
interaction g and the cavity-cavity inhomogeneities J\ 7^ J2 
constrain the excitation transfer, thus inhibiting a complete 
population tunneling. On the contrary, in the USC regime a 
much richer scenario appears. A complete photon transfer is 
generally allowed, even for finite values of the cavity-qubit 
coupling strength and for inhomogeneous hoppings, although 
there exists a specific regime for which the tunneling rate be- 
comes negligible. The complexity of the USC dynamics, that 
is generated by counter-rotating terms, manifests itself in the 
highly irregular time pattern of the observables, exhibiting a 
fractal behavior. Despite this fact, the physics beyond the 
RWA plays an important role for the enhancement of quan- 
tum state transfer. Finally, in the degenerate qubit case, the 
excitation transfer is regular and its period does not depend 
on the cavity-qubit coupling strength. 

The proposed scheme can be implemented with state-of- 
the-art current superconducting circuit technology or, alterna- 
tively, its quantum simulation may be considered. Further- 
more, it can be used as building block for realizing large 
controlled networks for quantum simulation and distributed 
quantum information processing involving the quantum Rabi 
model in all coupling regimes. 
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SUPPLEMENTARY MATERIAL 
State transfer: linear superposition 

In this section, we consider the state transfer of a Fock 
state linear superposition |0) and |1), containing zero and 
one excitations of the cavity mode, respectively. For sim- 
plicity, we have considered homogeneous cavity-cavity cou- 
plings, Jx = J-2, an d resonance condition, ui q = ui. The 
system is initialized in state \ipo) = \<fio) |0) |0) |g), where 
l^o) = p|0) + e y/l — p 2 |1), with p randomly chosen in 
[0,1]. We define the state transfer fidelity as F = Tr [p p(t)], 
which measures the overlap of the rightmost cavity state p(t) 
with the leftmost cavity initial state po = \<j>q) (4>o\. 

The results for different system parameters are shown in 
Fig. [6] On one hand, in Fig.[6h we have chosen a cavity-qubit 
coupling in the SC regime, g/u = 0.01, and two values for the 
cavity-cavity coupling J/ui = 0.01,0.1. Firstly, when both 
quantities are in the SC regime, the cavity-qubit interaction 
sets an upper bound for the maximum of the state transfer fi- 
delity (black line). Secondly, if the cavity-cavity coupling is in 
the USC regime, J/u) = 0.1, then almost complete state trans- 
fer occurs (red line). On the other hand, in Fig. |6b, we plot 
the fidelity F for a cavity-qubit coupling in the USC regime, 
g/bj = 0.2, and J/w = 0.01,0.1 for the cavity-cavity cou- 
pling. When g is one order of magnitude larger than J (black 
line), the transfer is completely inhibited. Contrariwise, there 
is an enhancement in the state transfer when the cavity-cavity 
coupling is in the USC regime (red line). 
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FIG. 7. (Color online) Transfer of a coherent state of amplitude a = 
1 along a three-cavity array (g = 0). (a) Average photon number 
for the leftmost (blue), central (green) and rightmost (red) cavities, 
(b) State transfer fidelity over time. In the case of qubit absence, the 
coherent state is fully transferred. Observe that the fidelity at time 
zero has a finite value F ~ 0.38, which corresponds to the overlap 
between the vacuum state and the considered coherent state. 



State transfer: coherent state 

In this section, we report the transfer properties of coher- 
ent states. We numerically simulated the system dynamics 
starting from the initial state |V>o) = |0o) |0) |0) |g), where 
\4>q) = \a) is a coherent state of amplitude |a|. For the sake 
of clarity, we first show how a coherent state would transfer 
along a three-cavity arrays in absence of the qubit interaction. 
This corresponds to setting g = in our model. Such dy- 
namics is shown in Fig.lTlfor a coherent state with a = 1. In 
FigfTL we show the average photon number for the three cav- 
ities, while Fig. fjh contains the state transfer fidelity. These 
figures show that the coherent state crosses the central cavity 
and it is recomposed in the rightmost one. The fast-oscillating 
behavior of the fidelity is due to relative phase rotation of dif- 
ferent coherent state components. We observe also that the fi- 
delity oscillates around a value slightly smaller than F = 0.4 
due to the finite overlap between the coherent state \a) and the 
vacuum state |0). 

Let us now consider our full model, composed of three cav- 
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FIG. 6. (Color online) (a) State transfer fidelity F = Tr [pop(t)] 
over time for different system parameters. We define \<f>o) = p |0) + 
e i8 y/l-p 2 \l}, with p = 0.2 and 6 = 0.63. In the SC regime, 
that is when g — J = 0.01 (black line), the state transfer fidelity is 
bounded: it cannot reach 1 as far as g 7^ 0. (b) When counter-rotating 
terms are involved in the dynamics (J > 0.1 and/or g > 0.1), the 
fidelity can be close to 1 even when J < g (red line). 



FIG. 8. (Color online) (a) State transfer fidelity over evolution time, 
for a coherent state of amplitude ex, — 1. The system parameters are 
given by g = 0.02 and J = 0.01 (SC regime), (b) State transfer 
fidelity over evolution time, for a coherent state of amplitude a — 1. 
The system parameters are given by g = 0.2 and J = 0.1 (USC 
regime). The maximum value that the fidelity can reach is not di- 
rectly bounded by the cavity-qubit interaction. 



ities and a qubit interacting with the central one. The trans- 
fer of coherent states follows the same general rules reported 
for the case of the linear superposition state, even if the time- 
dependent fidelity has a fast-oscillating behavior. In the SC 
regime, the cavity-qubit interaction limits the maximum value 
that the fidelity can reach. FigurelSh shows the time evolution 
of the state transfer fidelity in the case in which J = 0.01 and 
g = 0.01, where complete transfer is not allowed. In con- 
trast, in the USC regime, high values of the fidelity can be 



reached also when g is larger than J. Figure IHb shows the 
coherent state transfer dynamics, for a case in which both the 
cavity-cavity and the cavity-qubit interactions are in the USC 
regime. The plot shows the fidelity F over evolution time for 
a coherent state of amplitude |a| = 1. In this case, the fidelity 
can reach F = 0.9, also if g = 2 J, 

Notice that the considered cases, linear superpositions and 
coherent states, show the main features of the state transfer 
when considering the physics beyond the RWA. 



